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THEORETICAL  PREDICTIONS  FOR  SCATTERING  FROM  ELASTIC  SOUD  CYLINDERS  WITH 
HEMISPHERICAL  CAPS  AND  COMPARISON  WITH  DATA 

M.F.  WERBY  and  C.E.  DEAN 

Naval  Oceanographic  and  Atmospheric  Research  Laboratory,  Numerical  Modeling  Division,  Building 
1100,  Siennis  Space  Center,  MS  39529-5004,  USA 


The  extended  boundary  condition  equations  due  to  Waterman  have  led  to  a  T-matiix  that  can  describe 
scattering  6om  axially  symmetric  targets  over  a  broad  frequency  range.  Recent  interest  in  scattering 
from  cylinders  with  hemispherical  end  caps  has  prompted  us  to  analyze  several  examples  from  Le 
Havre.  All  phenomena  pre^cted  by  the  T-matrix  calculations  are  confirmed  by  the  data.  In  addition, 
bending  resonances  are  predicted  by  the  theoretical  calculations  and  compared  with  pr^ctions  from  bar 
theory.  A  series  of  predictions  are  made  for  various  materials  and  aspect  ratios  to  aid  future 
experimental  work. 

Les  Equations  des  conditions  dtendues  aux  limites  de  Waterman  nous  conduit  &  une  matrice  de  iransfert 
qui  peut  ddcrire  la  diffusion  par  des  cibles  h.  symdtiie  axiale  pour  un  grand  spectre  de  fri6quences. 
L'intdrfit  recent  pour  la  diffusion  par  les  cylindres  aux  bouts  hihnisphdriques  nous  incite  ^  I'analyse  de 
plusieurs  exemples  provenant  du  Havre.  Tous  les  phdnombncs  piddits  par  les  calculs  de  matrices  de 
transfert  se  confirnKnt  par  les  donndes.  De  plus,  les  resonances  de  flexion  se  predisent  par  les  calculs 
thd^ques  et  se  comparent  aux  predictions  de  la  theorie  des  banes.  Une  sdrie  de  predictions  se  fait  pour 
des  matieres  diverses  et  pour  des  rapports  divers  de  longueur  k  largeur  de  la  cible  pour  aider  les  travaux 
experimentaux  futurs. 

1.  Introduction 

Back  scattered  echoes  yield  interesting  information  when  scattering  from  elastic  targets,  particularly  in 
the  resonance  region.  Generally  one  has  source  and  receiver  at  the  sanoe  location  for  a  fixed  orientation.  If 
one  adjusts  die  orientation  over  a  suitable  angular  region  then  it  is  possible  to  ascertain  certain  symmetries 
of  the  target  if  they  exist  particularly  if  one  varies  the  frequencies  over  a  suitable  range  of  resonances.  If 
one  observes  axial  syrnriKtry  then  it  is  possible  to  obtain  both  the  dimensions  of  the  object  and  the  a^iect 
ratio  of  the  object  (ratio  of  length  to  width).  Once  this  information  is  obtained  it  is  also  possible  to  obtain 
information  on  whether  the  target  is  a  solid  or  a  shell  and  in  some  cases  what  the  shell  thickness  is  as  well 
as  the  material  properties  of  the  target  This  is  assuming  that  the  target  is  in  a  “free”  environment;  that  is, 
the  boundaries  of  the  target  are  not  a  factor  in  calculation.  Tune  domain  responses  for  specific  pulse  types 
also  yield  information,  and  it  is  easy  to  see  how  a  series  of  questions  can  form  the  basis  of  a  scenario  ^t 
can  rule  out  certain  targets  or  lead  to  a  probability  (confidence  level)  that  specific  targets  are  present 

Some  interest  has  arisen  for  the  case  of  scattering  from  elastic  solid  cylinders  with  hemisphmcal  end- 
caps.  We  have  performed  a  considerable  number  of  calculations  for  sphmids  with  aspect  ratios  out  to  10 
for  six  materials.  In  this  work  we  will  focus  on  scattering  from  a  tungsten  carbide  (WQ  cylinder  with 
hemispherical  end-caps  and  an  aspect  ratio  of  2.  We  wUl  describe  c^culations  for  scattering  end-on, 
broadside,  and  at  30°,  45°,  and  60°  relative  to  the  axis  of  symmetry.  A  brief  discussion  of  the 
methodology  and  physical  notions  that  pertain  to  the  final  analysis  is  presented  in  the  next  section  with  the 
filial  section  containing  the  analysis. 
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2.  Formal  expressions 

The  theoretical  calculations  presented  here  are  performed  using  noimal  mode  theory  for  the  spherical 
cases  and  by  the  extended  boundary  condition  (EBC)  or  T-mathx  method  fen'  the  nonsphetical  cases.  The 
EBC  (T-Matrix)  methods  have  been  publish^  elsewhere  so  we  will  not  present  the  details  here.  The 
boundary  conditions  and  pertinent  points  are,  however,  indicated  below.  The  boundary  conditions  at  a 
fluid-elastic  interface  are 


n«U+  =  n«U_,  n*t+  =  n»t_,  nxL  =  0,  (1) 

where  n,  U,  and  t  are  the  outward  unit  normal  to  the  object,  the  displacement  vector,  and  the  traction  on 
the  surface  of  the  scatterer,  respectively.  The  (-)  sign  indicates  an  exterior  (interior)  point  relative  to  the 
elastic  object.  The  T-matrix  relates  the  known  partial  wave  components  (An)  of  the  incident  field  (here  a 
plane  wave)  with  the  unknown  partial  wave  components  (fn)  of  the  scattered  field  as  follows: 

^n  =  ^  ^nn  *  ^n  •  (2) 

The  T-matrix,  which  is  derived  from  the  extended  boundary  condition  method,  is  symmetric  and  for 
nonlossy  objects  obeys  the  generalized  optical  theorem.  These  factors  can  be  built  into  the  formalism  or  T 
can  be  checked  to  assure  that  all  the  prop^es  are  satisfied,  and  that  convergence  has  been  achieved.  The 
form  function  is  obtained  from  the  T-matrix  and  the  incident  field  from  the  expression: 

f«(e)  =  2(kl72)  II  [(2n  l)(2n’  -t- 1 )]  1/2  jn  -  n'  P„.(cos  6)  An  I .  (3) 

The  rationale  behind  the  standing-wave  argument  for  determining  resonances  is  based  on  die  fact  that  for 
special  cases  one  can  equate  the  measure  P  around  the  perimeter  of  the  object  (in  the  direction  of  the 
circumnavigating  field)  with  a  half-integral  number  of  wavelengths.  For  the  case  (here)  for  elastic  solids 
we  are  close  to  a  rigid  background  and  one  has 

P  =  Xp(N-^l/2),  (4) 

where  N  =  2,  3,  4, ....  By  conservation  of  frequency  one  can  show  that 

Xp  =  (Vp/Vw)  •  (5) 

Since  kU2  =  2]tLAw  •  tbis  leads  to  the  standing  wave  condition  for  a  resonance 

[kLQlRes  =  L7i(Vj/V^)(N  +  l/2)/P .  (6) 

Uberall,  et  al.  have  shown  that  broadside  resonances  for  cylinders  can  be  obtained  from  the  sinqile 
expression 


rkU2]Res  =  (UD)(VpA^w)(N  +  1/2).  (7) 

Here,  we  may  use  the  phase  velocity  for  the  spherical  case  in  place  of  Vp.  Thus  with  simple  expressions 
one  can  predict  resonance  locations  for  both  end-on  and  broadside  resonances.  At  oblique  incidence  one 
can  not  only  excite  surface  resonances  but  also  bending  resonances.  It  has  been  shown  that  those 
resonances  are  related  to  bending  or  flexural  modes  which  are  also  predictable  using  Timoshenko  theory 
for  a  finite  circular  cylinder  (which  sqiproximates  a  spbertxd).  The  expression  for  the  resonance  location  in 
that  case  is  as  follows: 

fkL/2]Res  =  7t2/8V2  (Cl/Cw)(2n  -  1  )2/(a2  (a2  -k  3r7t2(n  -  1  )2/8)  1/2] ,  (8) 

where  a  designates  the  aspect  ratio  and  F  is  a  material  parameter. 
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3.  Analysis  of  results 

We  briefly  analyse  the  results  for  scattering  from  a  tungsten  carbide  cylinder  with  hemispherical  end- 
caps  for  a  range  of  dimensionless  half  length  size  parameter  kL/2  =  0  to  20  where  L  is  the  length  of  the 
object  and  k  the  wave  number  of  the  incident  field.  Often  one  chooses  r,  the  radius  of  the  target,  in  place 
of  the  half  length  L/2;  in  that  case  kr  would  range  from  0  to  10.  We  examine  results,  end-on,  broadside, 
and  at  oblique  angles  of  incidence.  In  the  future,  we  will  examine  results  out  to  klV2  =  30,  and  we  will 
subtract  a  rigidbackground  from  the  form  funcdon  leaving  only  the  resonance  response.  It  is  to  be 
remembered  that  we  expect  to  observe  Rayleigh  and  whispering  gallery  resonances  for  both  end-on  and 
broadside  incidence. 

Fig.  1  illustrates  the  backscattered  case  for  end-on  incidence.  One  observes  the  familiar  Franz  or 
creeping  wave  pattern  out  to  about  kL/2  =8  where  a  strong  (2, 1)  Rayleigh  resonance  is  present.  One  also 
observes  the  (3,1)  and  (4,1)  resonances  resonances  at  about  klj2  =  14  and  19,  respectively.  Whispering 
gallery  resonances  are  observed  at  intermediate  values.  Scattering  at  oblique  angles  relative  to  the  axis  of 
symmetry  yields  a  more  complicated  pattern.  This  is  to  be  expected  because  at  oblique  angles  resonances 
can  be  excited  ctxresponding  to  standing  waves  circumnavigating  the  largest  geodesic  (a  meridian)  and  the 
shortest  (parallel  to  the  equator  of  the  target).  Further,  flexural  or  bending  modes  can  also  be  excited  at 
oblique  angles.  Fig.  2  illustrates  the  case  for  backscattering  with  broadside  incidence.  There  we  see  a 
weak  Franz  wave  effect  characteristic  of  cylinders  as  well  as  a  bending  mode  at  about  kL/2  =  5.5  and  the 
end-on  Rayleigh  resonance  at  about  kL/2  =  8.  If  we  analyze  these  results  by  analogy  with  backscattering 
from  spheroids  then  the  first  major  dip  at  about  kL/2  =  13.5  is  the  first  broadside  Rayleigh  resonance 
while  the  second  dip  just  to  the  right  of  the  (2,  1)  resonance  is  the  (2,  2)  end-on  whispering  gall^ 
resonance  followed  by  two  broadside  whispering  gallery  resonances  and  finally  the  (3,  1 )  bre^side 
resonance  at  about  kL/2  =  20.  Fig.  3  has  the  same  features  as  the  broadside  scattering  case  including  the 
end-on  resonances  but  with  an  additional  effect  caused  by  bending  modes.  The  lowest  bending  mode  (n  = 
2)  occurrs  at  about  kL/2  =  2.75  and  the  second  mode  (n  =  3)  at  about  \dJ2=  5.5.  They  are  manifested  as 
sharp  peaks  and  are  less  pronounced  for  n  >  3.  It  has  been  shown  that  these  resonances  are  predicted 
using  Timoshenko  beam  Aoory.  Figs.  4  and  5  correspond  to  scanering  at  30°  and  60°  relative  to  the  axis 
of  symmetry.  All  features  are  present  in  these  cases  with  the  bending  mode  being  strongest  for  the  60° 
case. 

4.  Summary 

These  calculations  are  part  of  on-going  research  on  scattering  from  axially  symmetric  targets.  Extended 
calculations  are  to  follow  in  whichthe  background  will  be  subtracted  to  manifest  the  resonances  more 
clearly.  It  is  felt  that  these  results  coincide  with  recently  published  experimental  results  from  Le  Havre. 
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Fig.  1.  Backscanoed  fomi  function  vs.  kL/2  for  end-cn  incidence. 
(R^ponse  en  Cr^uence  de  rtoodiffusan  pour  incidence  selon  I'axe  de  syrndtrie.] 


kL/2 

Fig.  2.  Form  function  vs.  kL/2  for  broadside  incidence. 

(Rdponse  en  fodquence  de  idtrodiffiision  poor  incidence  i  90°  de  I'axe  de  symtoie.] 


Fig.  3.  Form  function  vs.  kL/2  for  incidence  a  45°  to  the  axis  of  symmeirie. 
[Rdponse  en  fidqoence  de  rdtrodiffiision  pour  incidence  k  45°  de  I'axe  de  symdoie.) 


